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ABSTRACT 

We present deep J and K s band photometry of 20 high redshift galaxy clusters between z = 0.8 — 1.5, 
19 of which are observed with the MOIRCS instrument on the Subaru Telescope. By using near- 
infrared light as a proxy for stellar mass we find the surprising result that the average stellar mass 
of Brightest Cluster Galaxies (BCGs) has remained constant at ~ 9 x lO n M0 since z ~ 1.5. We 
investigate the effect on this result of differing star formation histories generated by three well known 
and independent stellar population codes and find it to be robust for reasonable, physically motivated 
choices of age and metallicity. By performing Monte Carlo simulations we find that the result is 
unaffected by any correlation between BCG mass and cluster mass in either the observed or model 
clusters. The large stellar masses imply that the assemblage of these galaxies took place at the same 
time as the initial burst of star formation. This result leads us to conclude that dry merging has 
had little effect on the average stellar mass of BCGs over the last 9 — 10 Gyr in stark contrast to 
the predictions of semi-analytic models, based on the hierarchical merging of dark matter haloes, 
which predict a more protracted mass build up over a Hubble time. We discuss however that there is 
potential for reconciliation between observation and theory if there is a significant growth of material 
in the intracluster light over the same period. 
Subject headings: galaxies: clusters — galaxies: evolution — galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 

Brightest Cluster Galaxies (BCGs) are the most lumi- 
nous objects in the Universe in terms of stellar light and 
appear to be a separate population from bright ellipticals 
(Sandage 1972, 1976; Bhavsar & Barrow 1985; Oegerle 
& Hoessel 1991; Postman & Lauer 1995; Bernstein & 
Bhavsar 2001; Bernardi et al. 2007; von der Linden et al. 
2007; Vale & Ostriker 2008; Lin et al. 2009). They re- 
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side in the deep potential wells of the cores of rich galaxy 
clusters, thought to descend from the first regions where 
mass began to accumulate after the Big Bang. Their 
luminosities and unique environments make them ideal 
candidates for the testing of the mass build up in galaxies 
across a large fraction of the Hubble time. 

A number of studies have attempted to constrain the 
formation epoch and evolution of BCGs by compar- 
ing their K band Hubble diagram to a range of stellar 
population models (e.g. Aragon-Salamanca et al. 1998; 
Collins & Mann 1998; Nelson et al. 2002). BCGs have 
been shown to follow passive evolution out to moderate 
redshifts but then differing results are seen at z > 0.5, 
with some groups seeing a continuation of the passive 
trend (e.g. Collins & Mann 1998) while others claim the 
high redshift BCGs are fainter and therefore less massive 
than their local counterparts (e.g. Aragon-Salamanca 
et al. 1998). This is explained by a dependence of BCG 
luminosity on the mass of its host cluster (Burke et al. 
2000), with most studies now agreeing that the evolution 
of BCGs in massive clusters can be described as passive 
to z ~ 0.8 (Stott et al. 2008; Whiley et al. 2008). How- 
ever there is still some debate over the recent merging 
activity of BCGs particularly in lower mass clusters and 
Brightest Group Galaxies with a number of studies iden- 
tifying major merger candidates (Mulchaey et al. 2006; 
Rines et al. 2007; Tran et al. 2008). 

The favoured model for galaxy formation and evolu- 
tion is via the hierarchical merging of dark matter haloes 
(e.g. Davis et al. 1985). In this model the mass of 
a galaxy gradually increases as it merges with neigh- 
bouring systems. A major development in the field has 
been the advent of large cosmological simulations such as 
the Millennium N-body Simulation (Springel et al. 2005) 
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which models this hierarchical mass build up of dark 
matter haloes in a co-moving box 500/i _1 Mpc (where 

h = iVlOO kms^Mpc" 1 ) on the side. 

Semi-analytic models are commonly used to describe 
the complex baryonic physics of galaxies in the context of 
the merger histories of dark matter halos within N-body 
simulations (e.g. Bower et al. 2006; De Lucia & Blaizot 
2007). These models are an efficient way to describe the 
competing processes affecting baryonic matter such as 
those that trigger or suppress star formation i.e. gas rich 
galaxy merging and feedback from active galactic nuclei 
or supernovae. The output of observables, such as galaxy 
magnitudes, from the semi-analytic models has proved to 
be valuable for astronomers as the modelled systems can 
be compared directly with measurable quantities. For 
this reason semi-analytic models can be a useful tool 
for making mock catalogues, predicting the behaviour 
of galaxies, testing cosmological theory and assessing the 
feasibility of telescope observations. 

One important advantage of using BCGs to study 
galaxy evolution is that their theoretical counterparts 
can be easily and unambiguously identified as the cen- 
tral massive galaxies in mock clusters at the same red- 
shift as the real systems. De Lucia & Blaizot (2007) 
presented the evolution of BCGs over a Hubble time, in 
a semi-analytic model based on the Millennium Simula- 
tion. This paper predicted that the stellar population 
of BCGs forms early, with 50% of the stellar mass in 
place by z = 5 and 80% by z = 3. However this early 
star formation takes place in separate components that 
gradually assemble into the BCGs seen in the local Uni- 
verse, mainly through dry mergers that do not trigger 
additional star formation. So for example, at z = 1.5 the 
sum of the stellar mass in all sub-components is over 90% 
of the mass of the fully assembled BCG at z — 0, whereas 
the stellar mass in the main progenitor is on average only 
20% of the galaxy's final mass. 

From an observational standpoint, a significant num- 
ber of high redshift galaxy clusters {z > 1) have been dis- 
covered in recent years with X-ray surveys (e.g. Rosati 
et al. 2004; Mullis et al. 2005; Stanford et al. 2006; Bre- 
mer et al. 2006; Lamer et al. 2008). The XMM Clus- 
ter Survey (XCS, Romer et al. 2001; Sahlen et al. 2009) 
is one such project performing a serendipitous survey 
to discover clusters in the XMM-Ncwton archive. The 
main goals of the XCS are to constrain cosmological pa- 
rameters, measure the evolution of the hot gas through 
analysis of the cluster scaling relations and study galaxy 
evolution in clusters since the high mass cluster cores 
are thought to be environments comparable across all 
epochs. Furthermore X-ray luminosity and temperature 
are excellent proxies of cluster mass and enable us to 
directly compare the properties of real cluster galaxies 
with mock galaxies in similar cluster halo environments. 
With the advent of large wide-field optical and infrared 
surveys it is also possible to photometrically select high 
redshift galaxy clusters based on the properties of their 
constituent galaxies, further increasing the number of 
known z > 1 clusters (e.g. Barrientos et al. 2004; Swin- 
bank et al. 2007; Eisenhardt et al. 2008). 

In this paper we present near-infrared observations of 
20 high redshift galaxy clusters (0.8 < z < 1.5), de- 
rive stellar masses for their BCGs and compare to the 



latest semi-analytic models based on the Millennium N- 
body Simulation (De Lucia & Blaizot 2007; Springel et al. 
2005). We test how robust our results are to different star 
formation histories produced by three independent stel- 
lar population codes: Bruzual & Chariot (2003); Maras- 
ton (2005) and BaSTI (Pietrinferni et al. 2004; Percival 
et al. 2009). This work builds on the result of Collins 
et al. (2009) in which we found that the average stellar 
mass of the BCGs in 5 of the highest redshift galaxy clus- 
ters is not significantly different to that of BCGs in the 
local Universe. 

Unless otherwise stated a Lambda Cold Dark Matter 
(ACDM) cosmology (Q M =0.3, n Vac =0.7, H =70 km 
s _1 Mpc -1 ) is used throughout this work. 

2. CLUSTER SAMPLE 

Table 1 details our sample of 20 of the most distant, 
spectroscopically confirmed, galaxy clusters, including 
the highest redshift X-ray selected cluster XMMXCS 
J2215.9 - 1738 (Stanford et al. 2006; Hilton et al. 2007, 
2009). The sample consists of clusters discovered by var- 
ious X-ray surveys and several selected by optical meth- 
ods that show extended X-ray emission (see Table 1). 
The clusters all have spectroscopically confirmed red- 
shifts in the range 0.8 < z < 1.5 and have X-ray lu- 
minosities of 1 < Lx < 19 x 10 44 ergs _1 . 

To anchor our analysis to the local Universe we also 
include a low redshift comparison sample (0 < z < 0.3, 
0-7 < L x < 20.0 x lO^ergs- 1 ) published in Stott et al. 
(2008). This sample is a good low redshift comparison as 
the clusters cover a similar range in mass (average mass 
at z < 0.1 is 6.8(±1.5) x 10 14 M Q ) to the low redshift 
haloes of the Millennium Simulation (average mass at 
z = is 7.5(±3.5) x 10 14 M Q ). 

2.1. Cluster mass 

A number of authors have identified a weak correlation 
between BCG mass and their host cluster mass (Edge 
1991; Collins & Mann 1998; Burke et al. 2000; Brough 
et al. 2008; Stott et al. 2008; Whiley et al. 2008) which 
does not change significantly with redshift out to z ~ 0.8. 
Therefore in order to compare measured and predicted 
BCG masses in a meaningful way it is necessary that our 
cluster sample is well matched to the masses of simu- 
lated clusters in the Millennium Simulation with which 
we are comparing. The clusters of interest from the sim- 
ulation are the 125 most massive systems in the redshift 
snapshots z — 0.76, z — 1.08 and z — 1.5, selected for 
comparison with observations (De Lucia & Blaizot 2007). 
Halo masses M200 are measured at a radius (i?20o) inside 
which the average mass density is 200 times the critical 
density of the Universe. 

We calculate M200 and the associated uncertainties of 
our sample based on the observational results for the 
mass-temperature (M-Tx) relation, which is preferable 
to the L^-based determinations used in Collins et al. 
(2009) due to the putative presence of cluster cooling 
cores. The cluster X-ray temperatures used are listed in 
Table 1. We parameterize the M-Tx relation as 
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TABLE 1 

The cluster sample 



1. CL J0152. 7-1357 

2. XLSS ,1022303.0 - 043622 b 

3. XLSS ,1022400.5 - 032526 b 

4. RCS J0439 - 2904 

5. 2XMM ,1083026 + 524133 

6. RX ,10848.9 + 4452 c 

7. RDCS ,10910 + 5422 

8. CL J1008.7 +5342 

9. RX J1053.7 + 5735 (West) 

10. MS1054.4 - 0321 

11. CL J1226 -j- 3332 

12. RDCS ,11252.9 - 2927 

13. RDCS ,11317 + 2911 

14. WARPS ,11415.1 + 3612 

15. CL J1429.0 + 4241 

16. CL J1559.1 + 6353 

17. CL 16044-4304 

18. RCS J162009 + 2929.4 

19. XMMXCS J2215.9 - 1738 d 

20. XMMU J2235.3 - 2557 
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a The stellar mass errors quoted include the photometric error 
but not the uncertainty in the stellar population (sec §4.2 for a full 
treatment). 

b Bascd on XMM archival data analysed by XCS (Sahlcn et al. 
2009) 

c Archival photometry (Yamada et al. 2002). 

d Bascd on X-ray analysis to appear in Hilton et al. in prep 




Fig. 1. — 7 X 7arcsecond postage stamp K s band images of our 19 Subaru MOIRCS observed BCGs numbered as in Table 1. RX J0848.9 
I- 4452 is absent as this is an archival BCG (Yamada ct al. 2002) for which we do not possess Subaru MOIRCS imaging. 



the standard Hubble parameter at redshift z. The M 50 o 
masses are converted to M200 using the standard NFW- 
profilc prescription by Hu & Kravtsov (2003) with a halo 
concentration parameter c = 5. Wc include a a = 10% 
uncertainty on the normalization M* (similar to typical 
expected uncertainties in mass estimation, e.g. Nagai 
et al. 2007) along with the estimated measurement un- 
certainties on the temperatures, and from these derive 
the uncertainties on the M200 values by a simple Monte 



Carlo simulation. 

We use the parameter values based on the Maughan 
(2007) derived M-Tx relation, using the "centre ex- 
cluded" estimates in their Table 1: 



login (M* 



-0.57 , a = 1.72 , /? = -0.82 , a^ T = 0.06 . 



This normalization agrees with the relatively well- 
established value M 500 (z = 0.05, T x = 5keV) w 3 x 
lO^/i^M©, e.g. Pierpaoli et al. (2001); Reiprich & 
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Bohringer (2002); Viana et al. (2003); Vikhlinin et al. 
(2006). The scatter is the same as in Mantz et al. (2009), 
which is consistent with Maughan (2007) although the 
latter does not derive the M-Tx scatter explicitly. 

Finally we note that cluster mass estimates based on 
the somewhat steeper M-Tx relation by Mantz et al. 
(2009) give similar results. 

Crucially these cluster masses cover the range of mas- 
sive haloes seen in the equivalent redshift snapshots of 
the Millennium Simulation. The De Lucia & Blaizot 
(2007) simulated cluster samples at z — 0.76, z = 1.08 
and z — 1.5 have mass ranges at these redshifts of 2.4 — 
13.6x 10 14 M Q , 1.5-9.8 x 10 14 M Q and 1.0-7.5 x 10 14 M 
respectively. The average mass of the combined high 
redshift simulated haloes is 2.6 (±0.1) x 10 14 M Q , com- 
pared to the average mass for our sample of 4.5 (±0.7) x 
10 14 M©. However, based on the known trend of BCG 
mass with M200 (see §4.1) this 60% difference in average 
cluster mass equates at most to a shift in the average 
BCG mass of around 10%, which is within the measure- 
ment uncertainties. The consequences of this mismatch 
in cluster mass between the high redshift observations 
and the simulation is discussed in §4.3. 

3. OBSERVATIONS AND DATA REDUCTION 

The observations were taken with the MOIRCS cam- 
era (Ichikawa et al. 2006) on the 8.2m Subaru telescope 
which provides imaging and low-resolution spectroscopy 
over a total field-of-view of 4' x 7' with a pixel scale 
of 0.117" per pixel. This is achieved by dividing the 
Cassegrain focal plane and then re-focussing the light 
through identical optics onto two HAWAII-2 2048 x 2048 
CCDs, each covering 4' x 3.5'. Observations were taken 
in 0.5" seeing on the nights of August 8th and 9th 2007 
and in 0.3" - 0.6" seeing on the nights of December 16th 
2008 and April 18th 2009 with the clusters centered on 
Detector 2. A circular 11-point dither pattern of radius 
25" was used for both bands to ensure good sky sub- 
traction. The modal integration times were 25 mins at 
J and 21 mins at K s although we observed some of the 
higher redshift clusters for 50% longer when scheduling 
allowed. These exposures reach a 5cr limiting magnitude 
of at least J ~ 22.5 and K a ~ 21.5 (Vega). 

The data were reduced using the external IRAF pack- 
age MCSRED. The data were flat fielded, sky subtracted, 
corrected for distortion caused by the camera optical 
design, and registered to a common pixel coordinate 
system. The final reduced images on which we per- 
formed the photometry were made by taking the 3<r (s.d.) 
clipped mean of the dither frames. The BCG photometry 
was extracted in an identical manner to the low redshift 
comparison sample from Stott et al. (2008) using SEx- 
tractor (version 2.5) MAG_AUTO magnitude, which is 
found to be within ~0.1 magnitude of the total for ex- 
tended objects (Martini 2001). As a test of this method 
we also performed large aperture (50 kpc) photometry on 
the BCGs finding the values to be consistent with those of 
MAG_AUTO to within 0.05 mag. This ability to exclude 
light from close neighbors by using MAG_AUTO ensures 
that we are not including additional flux that would bias 
our stellar masses to higher values. We choose a global 
photometric background over a local one to control for 
any bias introduced by a low surface brightness halo or 
intracluster light in the vicinity of the BCG. To calculate 



the colors of the BCGs the images are first matched for 
seeing using the IRAF task PSFMATCH and then we 
run SExtractor in dual image mode so that the K s band 
detections extract the J band catalogue with identical 
positions and apertures to ensure accurate color deter- 
mination. This photometry is performed within fixed 8 
kpc circular apertures at the cluster redshift. 

The photometry was calibrated to the Vega system us- 
ing a combination of standard star observations and the 
2MASS and UKIDSS catalogues. The typical photomet- 
ric errors are 0.01 and 0.08 for the standard star and 
survey calibrated data respectively. 

Additional archival BCG photometry is included for 
the cluster RX J0848.9 + 4452 (Yamada et al. 2002) 
which is a total K s magnitude measured with a large 
aperture and is thus comparable to our own photometric 
analysis. 

3.1. BCG selection 

The BCG selection for a cluster is usually obvious from 
visual inspection of the images as for such rich clusters 
they are the prominent galaxy closest to the X-ray cen- 
troid often with a cD-like profile, however we chose to 
formalize this by studying the tip of the red sequence in 
the color magnitude relation. For each cluster we iden- 
tified the red sequence with J — K s color and selected 
the brightest galaxy from the K s -h&nd magnitudes of 
all the red sequence galaxies within a projected distance 
of 500 kpc from the cluster X-ray centroid as for ap- 
proximately 95% of clusters the BCG lies within this ra- 
dius (Lin & Mohr 2004). The only non-obvious case is 
J2215. 9-1738 where the object identified as the BCG is a 
spectroscopically confirmed member lying 300 kpc from 
the X-ray centroid which is only marginally brighter than 
several others in the cluster. A full discussion of this 
identification and the properties of other candidates is 
presented in Hilton et al. (2009). K s band images of our 
high redshift BCGs are shown in Fig. 1. 

4. STELLAR MASS 

4.1. Initial stellar mass calibration 

To compare the stellar mass build up in our observed 
BCGs with those of the semi-analytic model we need to 
derive stellar mass from the K s band luminosity of the 
galaxies. We first do this by calculating absolute K s band 
magnitudes for our galaxies by choosing a model appro- 
priate to the stellar population of the modelled BCGs. 
The De Lucia & Blaizot (2007) semi-analytic model pre- 
dicts a star formation history for the BCGs in which 50% 
of the stars have formed by z = 5 and 80% by z — 3, al- 
beit in separate components that have not yet coalesced 
to form the final mass of the galaxy. To model this evolu- 
tion a stellar population from Bruzual & Chariot (2003) 
with a Chabricr (2003) initial mass function (IMF) is 
implicit in the semi-analytic model. In Fig. 2 we investi- 
gate the stellar population of our observed BCG sample 
by plotting the J — K s color against redshift and com- 
paring to two Bruzual & Chariot (2003) simple stellar 
population (SSP) models with a Chabricr IMF and So- 
lar metallicity. A composite stellar population (CSP) 
with an exponentially decaying star formation rate with 
r=0.9Gyr which mimics the average stellar population 
present in the simulation (i.e. 50% of the stars have 
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formed by z = 5 and 80% by z = 3) is also plotted. 
From this plot we can see that although there is some 
scatter in the BCG near-infrared colors at z ~1 they are 
consistent with the CSP model. It should be stressed 
that this color evolution only gives information about 
the stellar population (e.g. age and metallicity) and not 
the mass of the system. 

With the initial assumption that our observed BCGs 
have similar stellar populations to those of De Lucia & 
Blaizot (2007) we calculate their stellar masses using the 
mass-to-light ratios from the CSP and normalise these 
masses and the model BCG masses from De Lucia & 
Blaizot (2007) at z = 0. We discuss the consequences 
of this assumption in §4.2 and compare a larger set of 
stellar population models appropriate to our BCGs. 

By using this technique we find stellar masses for our 
BCGs in the range 3.45 - 16.63 x lO n M with a biweight 
scale value of 8.52(±1.00) x 10 11 M Q which is consistent 
with the local BCG mass of 8.99(±0.82) x 10 n M Q . So 
when considering our entire sample with stellar masses 
derived directly from the CSP model we find that on 
average there has been no significant change in the stellar 
mass of BCGs out to at least z — 1.5. The corresponding 
stellar masses from the simulation at z = 0.76, 1.08, 1.5 
are respectively 3.84 ± 0.14,2.91 ± 0.10, 1.92 ± 0.07 in 
units lO n M0. The results of this analysis can be seen in 
Fig. 3 in which we plot stellar mass vs cluster mass for 
our BCG sample (black filled) and those of De Lucia & 
Blaizot (2007) at four different redshift snapshots (z=0, 
0.76, 1.0 and 1.5 corresponding to cyan, green, pink and 
red squares respectively). From this plot it is clear that 
the average mass of the observed BCGs is significantly 
higher than that of the model BCGs for clusters of similar 
mass and redshift. 

We note here that the choice of semi-analytic model 
does not affect our results as the 'Durham models' (i.e. 
Bower et al. 2006) also give a near identical discrepancy 
to that seen here (Collins et al. 2009). 

From Fig. 3 we see that the highest mass BCG is found 
in the highest mass cluster which may suggest a link 
between stellar mass and environment as seen by Edge 
(1991); Collins & Mann (1998); Brough et al. (2008); 
Stott et al. (2008); Whiley et al. (2008). However for our 
sample of 20 the correlation between BCG mass (Mbcg) 
and M200 is weak with a power law exponent of 0.42±0.12 
and a Spearman's rank analysis indicating this correla- 
tion is significant to the 99% level. Determination of 
this correlation in the literature from larger samples at 
lower redshift typically show a small dependency; for ex- 
ample, Whiley et al. (2008) find a power-law exponent 
of 0.12 ± 0.03, with a Spearman's rank correlation sig- 
nificant to greater than the 99.9% level. However, as 
mentioned in §2.1 it is clear from Fig. 3 that there is 
some mismatch in cluster mass between the high red- 
shift clusters and the corresponding simulated halos; we 
discuss the effect of this is §4.3. 

4.2. Dependence on stellar population models 

Rather than rely on the potentially naive stellar mass 
calibration used in §4.1 we look now at a range of stellar 
population models and codes with physically motivated 
input parameters to study their effects on the stellar mass 
evolution result. For completeness we utilize three lead- 
ing independent stellar population codes in this analysis 
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Fig. 2. — J — K s color vs redshift for our high redshift sample 
(red) with lex error bars and the low redshift analogue sample of 
Stott et al. (2008) (black). Two Bruzual & Chariot (2003) solar 
metallicity, simple stellar population models are included with for- 
mation redshifts Zj = 2 (dashed) and Zf = 5 (solid) as well as the 
model chosen to mimic the star formation histories of the De Lucia 
& Blaizot (2007) semi-analytic model which forms 50% of its stars 
by z = 5 and 80% by z = 3 (dotted). 



namely: Bruzual & Chariot (2003); Maraston (2005) and 
BaSTI (Pietrinferni et al. 2004; Percival et al. 2009), all 
of which are now widely used for extragalactic astronomy. 
In Collins et al. (2009) we investigated mass-to-light ra- 
tios extensively at z = 1.3 and found that our result 
held for the majority of combinations of stellar popu- 
lation synthesis code, age and metallicity. The notable 
exception to this was for young and sub-solar metallicity 
SSPs generated by the code of Maraston (2005) which, 
because of a strong emphasis on the Asymptotic Giant 
Branch phase of stellar evolution, gives young stellar pop- 
ulations (~ 1 Gyr) red colors degenerate with old age and 
high metallicity models. 

Until recently reliable metallicity determinations were 
available for only a few local BCGs. However Loubser 
et al. (2009) examined the stellar populations of 49 BCGs 
in the local Universe with high signal to noise spectra, 
concluding that on average they have at least twice-solar 
metallicity ([Z/H] = 0.31 ± 0.17) and enhanced a ele- 
ments ([E/Fe] = 0.41 ±0.09), suggesting an intense burst 
of star formation and subsequent quiescence. We use this 
information to rule out the low metallicity models and 
repeat the mass-to-light ratio analysis of Collins et al. 
(2009), concentrating on the variation in the age of the 
stellar component at twice-solar metallicity and includ- 
ing a enhancement (available only for BaSTI). Due to 
the differing metallicity sampling and IMFs available for 
the three stellar population codes, we use: Bruzual & 
Chariot (2003) with 2.5 Z and Chabrier IMF (Chabrier 
2003); Maraston (2005) with 2.2 Z Q and Kroupa IMF 
(Kroupa 2001); and BaSTI with 2.0 Z Q and Kroupa IMF. 
We calculate the mass-to-light ratios derived from these 
models for galaxies at the average redshift of our sample 
z = 1.0 and ages corresponding to formation redshifts 
2 < z f < 5 (2.6-4.7 Gyr at z = 1). The results of 
this analysis can be seen in Fig. 4, with the Bruzual & 
Chariot (2003) models giving the highest stellar masses 
followed by BaSTI (scaled solar then a enhanced) and 
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Fig. 3. — BCG stellar mass (Mq) vs cluster mass (Mq) for our sample (filled black circles z < 1 and filled black squares z > 1 to 
demonstrate that there is no redshift dependence). The error bars in the cluster mass include both the X-ray observation and full mass- 
temperature scaling errors, whereas the stellar mass errors include the photometric errors not the uncertainty in the stellar population see 
§4.2 and Fig. 4. The colored squares represent the BCGs from the model of De Lucia & Blaizot (2007) at the redshift snapshots 2 = 0, 
0.76, 1.0 and 1.5 (cyan, green, pink, red). The open circles show the z < 0.1 BCGs from the low z comparison sample of Stott et al. (2008) 



the Maraston (2005) models giving the lowest. The av- 
erage stellar mass predicted by the semi-analytic model 
at z = 1.0 is 2.91(±0.10) x lO n M . Therefore results 
of the mass-to-light ratio analysis have 3c significance 
or greater for all stellar population models considered 
here. This confirms that our result from §4.1 is robust 
to the influence of reasonable, physically motivated, star 
formation histories generated by independent synthesis 
codes. 

4.3. Dependence on cluster mass 

Due to the nature of detecting high redshift clusters 
in flux limited X-ray surveys, the clusters in our sample 
are relatively high mass systems. Because of this there is 
some mismatch between the average cluster mass of our 
sample and the average halo masses of the simulation 
(see §2.1). Given the relationship between cluster mass 
and BCG mass, previously discussed, this may lead to 
an unfair comparison. To account for this we perform a 
bootstrap simulation using the observations and the Mil- 
lennium Simulation haloes. The details of this simulation 
are as follows: we select a cluster at random from our ob- 
served sample, we then pick a model BCG from a halo of 



similar mass to the observation using a random normal 
selection with a sigma equivalent to the error on the X- 
ray inferred cluster mass as listed in Table 1 . To account 
for the discrete nature of the simulated redshifts we inter- 
polate the model BCG stellar mass to the redshift of the 
observed cluster. This is procedure is repeated 10,000 
times allowing for replacement. The resulting distribu- 
tion has an average stellar mass of 3.9 ± 0.2 x 10 n Af Q 
which is ~ 5er away from the average mass of our sam- 
ple, 8.52(±1.00) x 10 n M Q . Thus the significance of our 
result does not decrease when cluster selection effects are 
accounted for. 

5. SUMMARY AND DISCUSSION 

We have demonstrated that the average stellar mass 
of BCGs in the highest redshift X-ray clusters is dis- 
crepant with those from similar mass dark matter halos 
in semi-analytic models based on the Millennium Simu- 
lation. Instead of the gradual build up of mass through 
dry merging predicted by De Lucia & Blaizot (2007), 
our observations suggest that the stellar mass in these 
objects has remained unchanged over the last 9-10 Gyr 
requiring a more rapid build up of stellar mass before 
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Fig. 4. — Average stellar mass (Mq) vs formation redshift for 
our high redshift BCG sample when using mass-to-light ratios de- 
rived from three stellar population codes Bruzual & Chariot (2003), 
BaSTI (Pietrinferni et al. 2004; Percival et al. 2009) scaled solar, 
BaSTI a enhanced and Maraston (2005) (solid, short dashed, dot 
dashed and dotted lines respectively) with fixed twice-solar metal- 
licity, evaluated at the average redshift of our sample (z = 1.0). 
The upper long dashed line represents the average local BCG stel- 
lar mass 8.99±0.82 X 10 11 Mq while the lower long dashed line rep- 
resents the average stellar mass of the De Lucia & Blaizot (2007) 
semi-analytic model at z = 1 (2.91 ± 0.10 X 1O U M ). 

z = 1.5 some 4 — 5 Gyrs after the big bang. 

The timescale for the mass assemblage is similar to the 
age of the component stars (2 — 3 Gyrs), a situation that 
appears to resemble classical monolithic collapse (Eggen 
et al. 1962; Larson 1974) rather than hierarchical forma- 
tion. To form a galaxy of stellar mass 10 12 M© over 4 
Gyrs requires a mass deposition rate of about 250 M© 
yr _1 and an efficient mechanism to feed the gas into the 
inner regions of the halo where it can form stars. Un- 
fortunately the merging process becomes inefficient for 
massive galaxies because merger induced shocks lead to 
heating as opposed to radiative cooling of the gas (Bin- 
ney 2004). One suggestion is that the early assembly of 
massive galaxies at z > 2 is driven by narrow streams of 
dense cold gas which penetrate the shock-heated region 
greatly increasing the efficiency of the gas deposition and 
associated star formation (Birnboim & Dekel 2003; Keres 
et al. 2005; Dekel et al. 2009). 

Alternatively, a deficiency may lie in the semi-analytic 
treatment of the physical processes in the densest envi- 
ronments during early hierarchical assembly; this con- 
tention is supported by the fact that current predictions 
are moderately consistent with observations of the evolu- 
tion of luminous red galaxies (Wake et al. 2006; Almeida 
et al. 2008), whereas our results, which focus on the most 
massive subset of this population, the BCGs, differ much 
more from the model predictions. 

From a theoretical point of view this area of research 
is constantly evolving due to the challenges posed by ob- 
servation and the increased availability of powerful com- 
puters. A recent high resolution simulation of a single 
cluster (Ruszkowski & Springel 2009) predicts that 50% 
of the BCG final mass in a massive (10 15 M Q ) cluster has 
built up by hierarchical merging in the last 9 billion years; 
decreasing the discrepancy between our findings and the- 
ory. As this is for one cluster and not a full cosmological 



simulation on the scale of the Millennium Simulation it is 
difficult to know whether this improvement is due to low 
number statistics or a better approximation of N-body 
and stellar-dynamical effects. 

One consequence of our result is that if merging is im- 
portant at all since z ~ 1, the evolution of BCGs must 
be dominated by minor rather than major mergers, since 
the stellar mass appears unchanged since this time. Our 
observations are broadly consistent with the relatively 
low number of dry major mergers found at low redshift 
(Liu et al. 2009) and the model predictions of Khochfar 
& Silk (2009) which show that only 10 - 20% of galaxies 
more massive than 6.3 x 10 10 M© have experienced dry 
major mergers within their last Gyr at any given red- 
shift z < 1. Numerical simulations find that the scale 
sizes of early galaxies can grow from dry minor merging 
by a factor 2 — 3 since z = 1, e.g. Naab et al. (2009) 
and suggest that the cD-like haloes of BCGs are formed 
late, resulting in the relatively recent departure of BCGs 
from the Kormendy relation for ordinary elliptical galax- 
ies (Ruszkowski & Springel 2009). 

One possibility which might help reconcile observations 
and theory is the growth of stellar material constitut- 
ing the diffuse intracluster light within the cluster cores. 
From dissipationless simulations of dark matter haloes, 
Conroy et al. (2007b) find that while BCGs do not ap- 
pear to evolve strongly at z < 1, the intracluster light 
surrounding such galaxies is growing substantially, with 
up to ~ 85% of the stars in the intracluster medium of 
present day clusters deposited at z < 1 (see also Willman 
et al. (2004); Muranteet al. (2007); Conroy et al. (2007a); 
White et al. (2007); Henriques & Thomas (2010)). This 
inside-out growth is broadly consistent with the dry mi- 
nor merging scenario for the local elliptical galaxy pop- 
ulation (Bezanson et al. 2009; Hopkins et al. 2010) re- 
quired to explain their rapid size increase since z ~ 2 
(van Dokkum et al. 2009, 2010). 

Obscrvationally, recent results have confirmed the 
overall importance of intracluster light: Rudick et al. 
(2006) find that the intracluster light constitutes ~ 10% 
of the entire cluster stellar light; a result that appears 
to hold even for non-cD clusters (Feldmeier et al. 2004). 
From surface photometry out to 300 kpc of 24 clusters 
at z < 0.1, Gonzalez et al. (2005) demonstrate that the 
outer cD component of BCGs traces the cluster poten- 
tial and has ^10 times the total luminosity of the inner 
elliptical profile. These results suggest that further work 
on the growth of stellar light in the outer parts of BCGs 
is required to provide a consensus. 
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